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The long-time relaxation of the average conductance in a 
2D mesoscopic sample is studied within the method recently 
suggested by Muzykantskii and Khmelnitskii and based on 
a saddle-point approximation to the supermatrix cr-model. 
The obtained far asymptotics is in perfect agreement with 
the result of renormalization group treatment by Altshuler, 
Kravtsov and Lerner. 

PACS numbers: 71.55.Jv, 72.15.Lh, 05.40.+j 

In the recent paper fj| , Muzykantskii and Khmelnitskii 
(MK) considered the relaxation phenomena in disordered 
conductors in the framework of the supersymmetric cr- 
model approach. They suggested a nice idea that the 
long-time asymptotics of the conductance G(t) is gov- 
erned by a non-trivial saddle point of the cr-model. Their 
original goal was to reproduce in a more direct way the 
result of Altshuler, Kravtsov and Lerner (AKL) pj, who 
found the logarithmically normal (LN) "tail" in the time 
dispersion in two and (2 + e) dimensions. However, MK 
found in 2D a different, power-law decay for moderately 
large times. They put forward a hypothesis that the LN 
asymptotics could hold for longer times. Here I will show 
that this is indeed the case, and that this result can be 
obtained via the method developed by MK. 

Following MK, I consider a 2D disk-shaped sample of a 
radius R. I will consider the unitary symmetry (broken 
time reversal invariance) in course of the calculations. 
For systems of the orthogonal and symplectic symme- 
tries, the treatment is completely analogous, and I sim- 
ply present the corresponding results in the end of the 
paper. The problem can be described by the cr-model 
with the action @ 
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Here Q(r) is 4 x 4 supermatrix field, D is the diffusion 
constant, v the density of states, u> the frequency, Str 
denotes the supertrace, and A = diag(l, 1,-1,-1). The 
saddle point equation of MK reads: 
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where V n denotes the normal derivative. 

We can consider the two leads attached to the disk 
boundary to be of almost semicircular shape, with rel- 
atively narrow insulating intervals between them. Then 
we can approximate the boundary conditions by using 
eq.(||) for all the boundary, as it was done by MK. In 
fact, in view of the logarithmic dependence of the saddle 
point action on R (see below), the result should not de- 
pend to the leading aproximation on the specific shape 
of the sample and the leads attached. With the rota- 
tionally invariant form of the boundary condition, the 
minimal action corresponds to the function 8 depending 
on the radius r only. We get therefore the radial equation 



6" + 9'/r + 7 2 sinh0 = ; < r < R 
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(the prime denotes the derivative d/dr) with the bound- 
ary conditions: 



9(R) = , 
0'(O) = 
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The condition (0) follows from the requirment of analyt- 
icity of the field in the disk center. 

Assuming that characteristic values of 8 satisfy the 
condition S » 1 (wc will find below the corresponding 
restriction on the time t), one can replace sinh# by e e /2. 
Eq.(||) can be then easily integrated, and its general so- 
lution reads: 
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with two integration constants C\ and C 2 - To satisfy the 
boundary condition (Q), we have to choose C\ = 2. Fur- 
thermore, the above assumption 8(0) 3> 1 implies that 
4C2/7 2 ^ 1- Therefore, the second boundary condition 
(||) is satisfied if C 2 — (4/7i? 2 ) 2 , and the solution can be 
written in the form 



A L 8 + j 2 smli = 
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where 8(r) is the "non-compact angle" parametrizing 
the cr-model field Q(r), Al is the Laplace operator and 
7 2 = iuj/D. It should be supplemented by the boundary 
conditions on the boundary with leads 



Using now the self-consistency equation of MK, 

2n / dr r(cosh8 — 1) = t/nu , 
Jo 
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(3) one finds 7 2 = &ir 2 v/t. Finally, the action 
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is equal on the saddle point (Q) to 
S ~ &iz 2 vD ln(iA) 



where A = 1/(v-kR 2 ) is the mean level spacing. Eq.(|l] 
coincides exactly with the result of MK. This consider- 
ation is valid provided 0'(r) < l~ x on the saddle point 
solution, which is the condition of the applicability of the 
diffusion approximation (here / is the mean free path). In 
combination with the assumption 8(0) 3> 1 this means 
that l«fA« (R/l) 2 - 

Now I consider the ultra-long-time region, t ;g> 
A^ 1 (R/l) 2 . In order to support the applicability of the 
diffusion approximation, we should search for a function 
0(r) minimizing the action with an additional restriction 
0' < Since the derivative has a tendency to in- 

crease in the vicinity of r = 0, the restriction can be 
implemented via replacing the boundary conditions (Q) 
by #'(r*) = 0, where the parameter r* will be specified 
below. The solution reads now: 
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The function 9(r) is ment as being constant within the 
vicinity \r\ < of the disk center. The condition 
6' < yields ~ IC. It is important to note that 
the result does not depend on details of the cut-off pro- 
cedure. For example one gets the same results if one 
chooses the boundary condition in the form d'(r*) = 1/1. 
The crucial point is that the maximal derivative 0' should 
not exceed 1/7. The constant C is to be found from the 
self-consistency equation ([To]) which can be reduced to 
the following form: 
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Neglecting corrections of the ln(ln •) form, we find 
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The action mWj is then equal to 



S ~it 2 isD{C + 2) 2 ln(R/r*) 



For the orthogonal and symplectic ensembles, the 
saddle-point equation (|^) has the same form, with the 
only difference that the action (Jll|) is multiplied by the 
factor (3/2, where j3 = 1,2,4 for the orthogonal, uni- 
tary and symplectic symmetries respectively. Combining 
cqs.(|ll]) and (|l^), we get thus for the long-time asymp- 
totics of the average conductance G(t) ~ e~ s in all three 
symmetry cases: 



G(t)~ 
G(t) ~ exp • 
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where g = 2ixvD is the dimensionless conductance per 
square in 2D and r is the mean free time. 

The far asymptotical behavior (eq.([i8|)) is of the LN 
form and very similar to that found by AKL (see eq.(7.8) 
in Ref . ) . It differs only by the factor 1 / g in the argu- 
ment of In 2 . It is easy to see however that this difference 
disappears if one does the last step of the AKL calcu- 
lation with a better accuracy. Let us consider for this 
purpose the intermediate expression of AKL (Ref. [|, 
eq.(7.11)): 
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where u 



„ } n In in the weak localization region in 
2D, which we are considering. 



Evaluating the integral 
( |l9| ) by the saddle point method, we find 
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where we have kept only the leading term in the expo- 
nent. Eq.(p0[) is in exact agreement with eq.(|T8|) for 0=1 
(AKL assumed the orthogonal symmetry of the ensem- 
ble). Therefore, the supersymmetric treatment confirms 
the AKL result and also establishes the region of its va- 
lidity. It is instructive to represent the obtained results in 
terms of the superposition of simple relaxation processes 
with mesoscopically distributed relaxation times ts: 
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Then we have from eqs.(|l7j), ( |18[ ) for the distribution 
function P(t^): 
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where ~ R 2 /D is the time of diffusion through the 
sample. 

For completness, we list also the results for quasi-lD 
and 3D systems. For a quasi-lD sample (wire) of the 
length L (which is assumed to be much shorter than the 
localization length £ = 2/3-nvD) the asymptotics read 



-^^ln 2 (tAU , iA»l (23) 



G(t) ~ exp 



(for (3 = 2 this isjust eq.(16) of MK). It is interesting 
to note that eq. (E3) has essentially the same form as the 
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asymptotical formula for G(t) found by Altshuler and 
Prigodin [Q for the strictly ID sample with a length much 
exceeding the localization length: 



G(t)~exp^ln 2 (t/r) 



(24) 



If we replace in eq. (p4|) the ID localization length £ 
'21 by the quasi-lD localization length £ = 2f3irvD, we 
reproduce the asymptotics ( p5| ) (up to a normalization 
of t in the argument of In 2 , which does not affect the 
leading term in the exponent for t — * oo). This leads us to 
make the following two conclusions. Firstly, this confirms 
once more the general conjecture that the statistical 
properties of smooth envelopes of the wave functions in 
ID and quasi-lD samples are identical. Secondly, this 
shows that the asymptotical "tail" (|2^) in the metallic 
sample is indeed due to "quasi-localized" eigenstates, as 
has been conjectured (l],^|-|(|. 

In 3D, the analysis proceeds along the same line as for 
the ultra-long-time region in 2D. This is essentially what 
has been done by MK in their consideration of the 3D 
case. The result at t ^> {kfl) 2 tr> (where fc/ is the Fermi 
momentum) reads: 



G(t) ~ exp{-S*(£)} , S(t) ~ (k f l) 2 In 3 
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In contrast to the 2D case, the exact numerical coefficient 
in the exponent in eq.(p5|) cannot be found within the 
diffusion approximation. 

We note in conclusion, that the obtained long-time 
asymptotics of the average conductance have a very sim- 
ilar form to the asymptotical behavior of the distribution 
function P(p) of local density of states (LDOS) ||. In 
both cases, the result is of the LN form in quasi- ID and 
2D, and of a somewhat different (though very similar) 
exp{-(fc/Z) 2 ln 3 (-)} form in 3D. As in the case of LDOS 
distribution |^| , we have found a perfect agreement with 
the result of renormalization group (RG) treatment || 
in 2D. I believe this agreement between the RG and 
supersymmetric treatments of G(t) and P(p) to be of 
considerable conceptual importance. To make this point 
clear, I remind the reader that one of the first achiev- 
ments of the supersymmetry method as applied to dis- 
ordered electronic systems was the detailed study of the 
Anderson metal-insulator transition on the (effectively 
infinite-dimensional) Bethe lattice M. The found non- 
power-law critical behavior seemed at first sight to be in 
contradiction with the scaling hypothesis and with the 
results of RG treatment. Since the solution of the Bethe 
lattice problem was exact, this apparent contradiction 
questioned the validity of the scaling and RG approaches. 
These doubts were supported by the fact the solution in 
H heavily relied on the non-compact structure of the su- 
persymmetric cr-model manifold and was dominated by 
the large values 9 3> 1 of the "non-compact angle" 9. 



On the other hand, the RG consideration is just a re- 
summation of the perturbative expansion and does not 
distinguish between the compact and non-compact ver- 
sions of the a-model. 

However, we have been able to show recently Jl(J that 
the exotic critical behavior found in || is the property of 
infinite-dimensional models only and transits to a power- 
law one for a finite value d < oo of the space dimension, 
in qualitative agreement with predictions of the scaling 
and RG approaches. Results of || and of the present 
paper show a perfect quantitative agreement of super- 
symmetry and RG methods when applied to the problem 
of asymptotical behavior of various distributions in the 
ensemble of mesoscopic metallic samples in the weak lo- 
calization region. This provides strong support to other 
results obtained within the RG approach in the weak 
localization region and in the vicinity of the Anderson 
transition 0. On the other hand, we see that the su- 
persymmetry method is in many cases able to reproduce 
results of RG treatment in a more elegant way. Further- 
more, it is not restricted like RG to the spatial dimension 
d = 2 and can be successfully applied to quasi-lD and 
3D systems as well. Besides the study of conductivity 
relaxation G(t) and LDOS distribution P{p) discussed 
above, I would like to mention in this context the recent 
progress in understanding of the statistical properties of 
eigenf unctions |0,O,^2[. Seeing that the two approaches 
are in amazingly good agreement, we can (depending on 
the problem considered) use any of them or even combine 
them to complete our understanding of the properties of 
mesoscopic disordered systems. 

I am grateful to V.E.Kravtsov and D.E.Khmelnitskii 
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